i 


/9 J 5-/4  /  W 


s/ 

AFWAL-TR--S6-4016 


ESTIMATING  THE  STATISTICAL  PROPERTIES  OF  CRACK  GROWTH  FOR  SMALL  CRACKS 


P.  W.  Hovey 
J.  P.  Gallagher 
A.  P.  Berens 


University  of  Dayton  Research  Institute 
300  College  Park  Drive 
Dayton,  Ohio  45469 


December  1980 


TECHNICAL  REPORT  AFWAL-TR-80-4016 

Final  Report  for  Period  September  1979  -  September  1980 


Approved  for  public  release;  distribution  unlimited 


MATERIALS  LABORATORY 

AIR  FORCE  WRIGHT  AERONAUTICAL  LABORATORIES 

AIR  FORCE  SYSTEMS  COMMAND 

WRIGHT- PATTERSON  AIR  FORCE  BASE,  OHIO  45433 


20070924053 


NOTICE 


When  Government  drawings ,  specificationsr  or  other  data  are  used  for  any  purpose 
other  than  in  connection  with  a  definitely  related  Government  procurement  operation, 
the  United  States  Government  thereby  incurs  no  responsibility  nor  any  obligation 
whatsoever ;  and  the  fact  that  the  government  may  have  formulated,  furnished ,  or  in 
any  way  supplied  the  said  drawings,  specifications ,  or  other  data,  is  not  to  be  re¬ 
garded  by  implication  or  otherwise  as  in  any  manner  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  rights  or  permission  to  manufacture 
use,  or  sell  any  patented  invention  that  may  in  any  way  be  related  thereto. 

This  report  has  been  reviewed  by  the  Office  of  Public  Affairs  (ASD/PA)  and  is 
releasable  to  the  National  Technical  Information  Service  (NTIS) .  At  NTIS,  it  will 
be  available  to  the  general  public,  including  foreign  nations. 


WALTER  H.  RE  I  MANN 
Acting  Chief 
Metals  Behavior  Branch 
Metals  and  Ceramics  Division 


Project  Engineer 
Metals  Behavior  Branch 
Metals  and  Ceramics  Division 


"If  your  address  has  changed,  if  you  wish  to  be  removed  from  our  mailing  list,  or  if  the 
addressee  is  no  longer  employed  by  your  organization  please  notify  AFWAL/MLLN  W-P  APB, 

OH  45433  to  help  maintain  a  current  mailing  list". 

Copies  of  this  report  should  not  be  returned  unless  return  is  required  by  security  considerations, 
contractual  obligations,  or  notice  on  a  specific  document. 


SECURITY  cl  OSSIFICATION  OF  this  PAGE  (Wh*n  Data  Entered) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

1.'  REPORT  NUMBER  2.  GOVT  ACCESSION  NO. 

AFWAL-TR-81 -4016 

3.  RECIPIENT'S  CATALOG  NUMBER 

*.  TITLE  (and  Subtitle) 

ESTIMATING  THE  STATISTICAL  PROPERTIES 

OF  CRACK  GROWTH  FOR  SMALL  CRACKS 

5.  TYPE  OF  REPORT  A  PERIOD  COVERED 

Final  Report 

Sept.  1979  -  Sept.  1980 

6.  PERFORMING  ORC.  REPORT  NUMBER 

UDR-TR -80-123 

7.  AUTHORfjj 

P.  W.  Hovey,  J.  P.  Gallagher,  and  A.  P,  Berens 

e.  CONTRACT  OR  GRANT  NUMBERf*) 

F3361 5-78-C-51 84 

9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

University  of  Dayton 

Research  Institute 

Dayton,  Ohio  45469 

10.  PROGRAM  ELEMENT,  PROJECT,  TASK 

AREA  A  WORK  UNIT  NUMBERS 

2418  -  Prog.  No. 

03  -  Task  No. 

06  -  Work  Unit  No. 

M.  CONTROLLING  OFFICE  NAME  AMO  ADDRESS 

Materials  Laboratory  (AFWAL/MLLN) 

Air  Force  Wright  Aeronautical  Laboratories 
Wriqht-Patterson  Air  Force  Base,  Ohio  45433 

12.  REPORT  DATE 

March  1981 

13.  NUMBER  OF  PAGES 

26 

U.  MONITORING  AGENCY  NAME  a  ADDRESS^/  different  from  Controlling  Office) 

15.  SECURITY  CLASS,  (of  tbia  report ) 

Unclassified 

15a.  DECLASSIFICATION 'DOWNGRADING 

schedule 

16  DISTRIBUTION  STATEMENTfonh^^porfj 


Approved  for  public  release;  distribution  unlimited. 


J  '7.  DISTRIBUTION  STATEMENT  (of  the  aba  tract  entered  in  Block  JG,  ft  different  from  Report) 


IB.  supplementary  notes 


19  KEY  WORDS  (Continue  on  toverae  aide  it  neceeaary  and  identity  by  block  number) 


20.  ABSTRACT  f Confirm®  on  raver ee  aide  it  neceeaary  and  identify  by  block  number) 

A  series  of  statistical  studies  were  conducted  on  two  sets  of  fatigue  crack 
growth  data  which  were  generated  under  variable  amplitude  loading  from  cracks 
that  initiated  in  7475-T7351  Aluminum  down  in  the  bore  of  typical  aircraft 
quality  prepared  holes.  The  study  considered  fatigue  crack  growth  (life) 
behavior  and  fatigue  crack  growth  rate  behavior.  Generally,  the  study  showed 
that  small  crack  growth  behavior  was  similar  to  that  observed  for  larger  cracks 
and  that  a  stress  intensity  factor  parameter  could  be  used  to  describe  the  mean 
trend  of  the  fatigue  crack  growth  rate  behavior  via  a  power  law  relationship. 


DD  ,^3  1473  EDITION  OF  1  NOV  6*  IS  OBSOLETE 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  r»Ti®n  Data  Entered) 


FOREWORD 


This  study  was  conducted  by  the  authors  during  the  time 
frame  from  September  1979  to  September  1980  under  USAF  contract 
F33615-78-C-5184  for  the  Air  Force  Wright  Aeronautical  Labora¬ 
tories/the  Materials  Laboratory.  Dr.  Theodore  Nicholas  of  the 
Metals  and  Ceramics  Division  of  the  Materials  Laboratory  was 
the  project  monitor  for  this  study  of  the  application  of  non¬ 
linear  fracture  mechanics  (NLFM)  parameters  to  the  study  of 
fatigue  crack  growth.  One  of  the  purposes  of  considering  the 
small  crack  data  presented  in  this  report  was  to  determine  if 
NLFM  parameters  were  required  to  characterize  small  crack 
behavior.  For  the  data  sets  considered,  the  data  were  adequately 
described  by  the  linear  elastic  fracture  mechanics  approach. 
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SECTION  1 
INTRODUCTION 


This  report  investigates  some  of  the  statistical  charac¬ 
teristics  of  fatigue  crack  growth  (FCG }  behavior  for  small  cracks. 
The  characteristics  covered  include  FCG  and  fatigue  crack  growth 
rate  {FCGR)  variability,  crack  length  serial  effects,  and  the 

relationship  of  FCGR  to  stress  intensity  factor  (K  ) .  Crack 

max 

length  serial  effects  can  be  described  as  the  relationship  between 
observed  FCGR's  for  adjacent  increments.  Because  the  material  is 
expected  to  be  similar  in  adjacent  increments,  the  FCGR’s  are 
also  expected  to  be  similar.  Besides  presenting  these  descriptive 
statistics,  a  comparison  is  made  between  large  and  small  cracks. 

Various  aspects  of  the  statistics  of  FCG  and  FCGR  have 
been  studied  by  others.  Mukherjee  and  Burns1  studied  the  use 
of  regression  models  relating  FCGR  to  the  stress  intensity  factor 
and  discussed  the  errors  associated  with  these  methods.  The 
effects  of  different  differentiation  methods  and  increment  sizes 
on  FCGR  variability  for  constant  stress  intensity  factor  conditions 

■  ,  2  7 

was  investigated  by  Artley,  et  al  .  Shaw  and  LeMay  recently 
compared  various  methods  of  determining  intrinsic  material  scatter. 
An  extensive  set  of  FCG  data  was  collected  for  constant  amplitude 
stress  conditions  by  Virkler,  et  al.4  in  a  study  aimed  at 
evaluating  various  statistical  distributions  to  FCGR  at  different 
stress  intensity  factor  levels.  As  part  of  their  study  Virkler, 
et  al.  also  attempted  to  correlate  the  variability  in  FCGR 
behavior  with  the  variability  in  FCG  life  behavior.  The  Virkler, 
et  al.  FCG  data  set  was  also  used  as  the  basis  for  a  recent 

5 

study  by  Berens,  et  al.  which  modeled  FCG  lives  from  FCGR  data 
derived  from  the  original  FCG  data  set.  Simulations  were  used  in 
the  Berens  et  al.  study  to  validate  the  model  that  was  developed. 
Some  results  developed  in  these  earlier  studies  will  be  compared 
to  results  generated  in  this  study. 
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The  data  for  the  current  study  were  obtained  from  a  study 
which  was  designed  to  evaluate  the  influence  of  the  quality  of  a 
fastener  hole  on  structural  durability  of  7475-T7351  Aluminum. 

While  this  study  investigated  a  wide  range  of  loading  and  fastener 
conditions,  only  two  drilling  machines  were  used:  The  Quackenbush 
Model  No.  QDA-100  and  the  Winslow  Spacematic  Model  No.  SI.  Besides 
being  different  manufacturers,  these  represent  two  different  pro¬ 
cedures.  The  Quackenbush  process  is  a  two-stage  drill  and  ream 
procedure  while  the  Winslow  process  includes  only  drilling. 

The  FCG  data  selected  for  the  current  study  were  taken  from 
no-load-transfer  specimens  which  had  properly  drilled  holes  and 
which  were  tested  under  a  F-16  fighter  load  spectrum.  The  FCG 
data  sets  were  further  censored  to  include  only  those  specimens 
exhibiting  FCG  through  the  crack  length  interval  from  0.0025  to 
0.0375  inches.  While  the  exact  shape  of  each  crack  is 
unknown,  the  investigation  revealed  that  typically  the  cracks  were 
growing  in  the  bore  of  the  hole  in  a  semielliptical  shape. 
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SECTION  2 

VARIABILITY  OF  LIFE 


In  the  Reference  6  study,  FCG  data  were  collected  periodically 
with  a  constant  time  increment  frequency  of  400  flight  hours  (the 
length  of  the  repeating  load  block) .  Interpolation  was  used  to 
create  a  set  of  FCG  data  with  a  constant  crack  growth  increment 
(Aa)  of  0.005  inches  through  the  same  interval.  A  constant  Aa 
facilitates  the  investigation  of  the  effect  of  stress  intensity 
factor  on  FCGR  and  of  serial  effects  on  FCGR .  The  interpolated 
data  were  normalized  to  zero  life  at  a  crack  length  of  0.0025 
inches  and  since  the  crack  grew  to  0.0375  inches  this  gives  seven 
increments  to  calculate  FCGR's. 

The  two  sets  of  FCG  data  are  presented  in  Figures  1  and  2. 
Figure  1  contains  the  data  obtained  from  specimens  prepared  with 
the  Quackenbush  tooling,  likewise,  Figure  2  contains  the  data  for 
the  test  series  where  the  hole  was  prepared  using  the  Winslow 
tooling.  The  code  names  for  these  two  data  sets  are  QPF  and  WPF. 

The  P  stands  for  proper  drilling  and  the  F  for  fighter  spectrum. 

The  Q  and  W  have  the  obvious  meanings  of  Quackenbush  and  Winslow. 

In  the  WPF  group,  as  shown  by  Figure  2,  the  FCG  behavior 
from  two  specimens  stand  out  from  the  rest.  These  two  specimens 
are  coded  WPF  7  and  WPF  11  in  the  data  report.  The  values  of 
the  test  statistics  (r  ratios)  are  0.76  and  0.81  for  WPF  7  and 
WPF  11,  respectively;  both  of  these  values  lie  above  the  99.5 
percentile  of  the  distribution  of  r  ratios.  If  the  fatigue  lives 
for  WPF  7  and  WPF  11  came  from  the  same  distribution  as  the  rest 
of  the  WPF  specimens  fatigue  lives,  there  would  be  less  than  one 
half  of  one  percent  chance  that  values  for  the  r  ratio  as  high 
as  those  observed  would  occur.  This  indicates  that  something 
different  happened  in  the  testing  of  specimens  WPF  7  and  WPF  11. 
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FLIGHT  HOURS  (F)  IN  THOUSANDS  OF  HOURS 

Figure  1.  Fatigue  Crack  Growth  Behavior  Exhibited  by  Open  Hole  Specimens, 
Prepared  with  the  Quackenbush  Equipment  and  Subjected  to  an  F-16 
Fighter  Wing  Stress  History  (Referred  to  as  the  QPF  Data  Set) . 
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FLIGHT  HOURS  (F)  in  thousands  of  hours 

Figure  2.  Fatigue  Crack  Growth  Behavior  Exhibited  by  Open  Hole  Specimens 
Prepared  with  the  Winslow  Equipment  and  Subjected  to  an  F-16 
Fighter  Wing  Stress  History  (Referred  to  as  the  WPF  Data  Set) . 


Although  there  is  no  record  of  abnormalities  in  the  test 
conditions  for  these  two  specimens,  they  will  be  removed  from 
subsequent  analyses.  The  main  purpose  of  this  study  is  to 
collect  background  data  for  the  development  of  a  model  of  FCG 
lives  based  on  FCGR.  Initial  modeling  efforts  require  a 
simplicity  which  cannot  be  obtained  with  outliers  present. 
Henceforth,  in  this  report  the  term  WPF  set  will  refer  to  the  set 
with  WPF  7  and  WPF  11  excluded,  unless  otherwise  indicated. 

Summary  statistics  for  "lifetimes,"  in  flight  hours,  for 
the  crack  to  grow  from  0.0025  inches  to  0.0375  inches  are  listed 
in  Table  1.  Three  sets  of  statistics  are  presented:  one  for  the 
QPF  group,  one  for  the  WPF  group,  and  one  for  the  WPF  group  with 
the  two  outliers  mentioned  in  the  previous  paragraph.  The  mean, 
standard  deviation,  and  coefficient  of  variation  of  both  lifetimes 
and  log  lifetimes  are  given  in  the  table.  The  QPF  set  seems  to 
grow  faster  than  the  WPF  set,  particularly  when  the  two  outliers 
are  ignored . 

Typically,  a  lognormal  distribution  is  used  to  describe  the 
distribution  of  lifetimes  or  of  fatigue  life  ratios.  Life  ratios 
are  the  ratio  of  predicted  life  to  actual  life.  For  these  data, 
life  ratios  were  calculated  using  the  mean  lifetime  for  the  set 
(QPF  or  WPF)  as  the  predicted  life.  Figure  3  contains  plots  of 
the  observed  cumulative  distribution  functions  on  a  lognormal 
probability  scale.  As  can  be  seen  in  Figure  3,  all  of  the  data 
fall  very  close  to  a  straight  line  which  is  indicative  of  a  log¬ 
normal  distribution  for  life  ratios.  This  can  be  substantiated 
with  the  Shapiro  Wilk  W  test.  This  is  a  test  for  the  applicability 
of  the  normal  distribution  and  when  applied  to  log  values  can  be 
used  to  test  for  the  lognormal  distribution.  The  W  statistics 
for  the  QPF  set  and  the  WPF  set  without  the  two  outliers  are 
0.9741  and  0.9462,  respectively.  In  this  test,  low  values 
represent  non-normality  (or  non  lognormality ) .  The  values  of 
W  for  the  QPF  and  WPF  without  outliers  sets  are  between  the  50th 
and  90th  percentiles  of  the  W  distribution,  implying  that  the 
assumption  of  a  lognormal  distribution  is  reasonable. 
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SUMMARY  STATISTICS  FOR  LIFETIMES  FOR  CRACK  GROWTH  FROM 
0.0025  INCHES  TO  0.0375  INCHES 
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CUMULATIVE  PROBABILITY,  PERCENT 


.2  4  12  4 


Figure  3.  Fatigue  Life  Ratios  Cumulative  Distributions. 
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An  F  ratio  test  can  be  used  to  test  for  the  equality  of 
standard  deviations  of  log  lives.  The  standard  deviation  of  a 
sample  is  the  square  root  of  the  sample  variance.  It  is  well 
known  that  the  ratio  of  variances  of  samples  from  two  different 
normal  populations  with  the  same  variance  has  an  F  distribution. 

If  the  observed  F  ratio  is  near  one,  then  it  would  seem  likely 
that  the  variances  are  the  same.  If  it  is  much  larger  or  smaller 
than  one,  the  variances  are  probably  different.  The  F  ratio 
(QPF  over  WPF)  for  comparing  the  QPF  set  to  the  WPF  set  without 
the  outliers  was  4.63.  We  would  expect  the  F  ratio  to  be  less 
than  2.65  about  ninety-five  percent  of  the  time.  Thus,  we  feel 
confident  that  the  variances  are  different. 

The  mean  of  log  life  ratios  is  always  close  to  zero,  simply 
because  of  the  way  they  are  calculated.  Therefore,  it  is 
impossible  to  see  any  differences  in  mean  log  lives  in  Figure  3. 

To  investigate  this  possibility,  the  t  test  was  used.  The  value 
of  t  for  testing  for  a  difference  in  means  for  the  QPF  set  and 
the  censored  WPF  set  was  4.90.  This  is  well  above  the  99.5  per¬ 
centile  of  a  t  with  32  d.f.  distribution  (2.741).  In  other  words, 
there  is  less  than  a  one-half  of  one  percent  chance  of  getting  a 
value  of  4.9  or  larger  for  t  if  the  two  means  were  the  same. 
Therefore,  it  is  fairly  certain  the  mean  log  lifetimes  are  different. 
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SECTION  3 

VARIABILITY  IN  FCGR 

The  variability  in  FCGR  for  the  short  crack  data  analyzed 
showed  similar  characteristics  to  FCGR  variability  observed  in 
studies  involving  longer  cracks.  The  mean,  standard  deviation, 
and  coefficient  of  variation  of  FCGR  are  listed  as  a  function 
of  crack  length  in  Table  2.  The  FCGR  data  were  created  using  the 
secant  method  on  0.005  inch  increments.  The  statistics  are 
calculated  for  each  crack  growth  increment  within  each  data  set. 

The  crack  length  (a)  given  is  the  midpoint  of  the  increment. 

The  only  readily  apparent  trend  is  that  the  mean  FCGR  increases 
with  the  crack  length.  On  closer  observation,  the  coefficients 
of  variation  for  the  QPF  specimens  are  noted  to  be  generally 
larger  than  those  of  the  WPF  group,  and  the  mean  FCGR's  for  QPF 
specimens  are  noted  to  be  generally  higher  than  for  WPF  specimens 
in  each  increment  considered. 

In  both  cases,  the  coefficients  of  variation  are  within  the 

bounds  seen  by  other  studies.  In  particular,  a  study  by  Artley, 

2 

et  al.  showed  coefficients  of  variation  for  secant  method 

differentiation  procedures  ranging  from  about  35%  for  a  0.008 

inch  increment  to  15%  for  a  0.034  inch  increment  as  shown  by 

Figure  4.  Figure  4  also  provides  a  comparison  between  several 

studies  including  the  Artley,  et  al.  study,  a  study  by  Virkler, 

et  al.4,  the  fastener  hole  quality  study,  and  a  study  by  Shaw 

and  LeMay^ .  A  trend  in  coefficient  of  variation  as  a  function  of 

increment  size  is  also  illustrated  in  this  figure.  The  materials 

2 

employed  in  these  studies  were  7075-T6  Aluminum  (Artley,  et  al.  ) , 

4 

2024-T3  Aluminum  (Virkler,  et  al.  ),  and  4140  Steel  (Shaw  and 

3 

LeMay  ) .  It  is  believed  that  the  accuracy  of  the  crack  length 

measurements  is  likely  to  be  the  cause  for  the  differences  seen. 

The  higher  measurement  accuracy  associated  with  electron 

fractography  would  explain  the  fact  that  the  coefficient  of 

variation  in  the  fastener  hole  quality  study  is  lower  than  that 

2  . 

observed  by  Artley,  et  al.  ,  who  employed  light  microscopes  (10-30x) . 
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MEAN,  STANDARD  DEVIATION,  AND  COEFFICIENT  OF 
VARIATION  OF  FCGR  AS  A  FUNCTION  OF  CRACK  LENGTH 
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SECTION  4 

SERIAL  CORRELATION 


When  observations  are  taken  close  together  either  spacially 
or  in  time,  there  is  the  possibility  that  they  are  related.  This 
is  referred  to  as  serial  correlation  or  autocorrelation.  In  this 
study,  observations  of  FCGR  are  taken  in  successive  increments 
of  crack  growth,  therefore  serial  correlation  is  suspected. 

Cross  correlations  are  used  to  investigate  serial  correla¬ 
tions  in  this  study.  Cross  correlation  is  a  term  used  to  refer 
to  the  correlation  coefficients  of  all  possible  pairs  of  variables 
from  a  specified  group  of  variables.  In  this  case,  the  seven 
variables  are  the  FCGR  values  from  the  seven  increments.  Cross 
correlations  are  presented  in  matrix  form  with  the  row  and  column 
of  an  entry  designating  the  two  variables  for  which  it  is  a  measure 
of  association.  Stated  differently,  the  correlation  coefficient 
for  the  second  and  fourth  variables  is  entered  in  the  fourth 
column  of  the  second  row.  If  the  full  matrix  were  presented, 
this  correlation  would  also  be  entered  in  the  second  column  of 
the  fourth  row.  However,  to  avoid  this  duplication,  usually 
only  the  upper  right-hand  corner  of  the  matrix  of  correlation 
coefficients  is  presented.  Tables  3  and  4  contain  the  correlation 
matrices  for  the  QPF  and  WPF  without  outliers  data  sets, 
respectively. 

Values  of  the  correlation  coefficient  range  between  -1  and 
1.  A  value  of  1  or  -1  is  considered  perfect  correlation  and 
indicates  that  one  variable  can  be  predicted  exactly  from  the 
other.  That  is,  if  the  FCGR  is  one  unit  above  the  mean  in  the 
first  increment,  the  FCGR  for  that  specimen  in  the  second  incre¬ 
ment  will  also  be  1  unit  above  the  mean.  A  value  of  0  for  the 
correlation  coefficient  indicates  no  relationship  between  the  two 
variables;  and  therefore,  the  FCGR's  in  the  two  increments  are 
independent.  Values  between  0  and  1  or  0  and  -1,  indicate  vary¬ 
ing  degrees  of  the  strength  of  the  relationship  between  the  two 
variables. 
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TABLE  3 

CROSS  CORRELATIONS  FOR  THE  QPF  DATA 
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CROSS  CORRELATIONS  FOR  THE  WPF  DATA  WITHOUT  OUTLIERS 
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Usually  when  the  variables  are  ordered,  as  these  are,  the 
correlation  coefficient  is  assumed  to  be  a  function  of  the  dis¬ 
tance  between  the  two  variables.  For  example,  it  would  be  assumed 
that  the  first  and  second,  the  second  and  third,  the  third  and 
fourth,  .  .  . ,  sixth  and  seventh  variables  all  have  the  same 
correlation  coefficient  since  they  are  all  1  step  apart.  Thio 
relative  separation  is  referred  to  as  the  lag  so  that  the  connected 
pairs  represent  a  lag  of  1.  For  a  lag  of  2,  the  pairings  would 
be:  first  and  third,  second  and  fourth,  third  and  fifth,  fourth 

and  sixth,  and  fifth  and  seventh. 

In  Tables  3  and  4  all  the  lag  1  pairings  are  in  the  second 
diagonal  of  the  matrix.  The  first  diagonal  contains  all  l's 
since  it  represents  the  correlation  of  each  variable  with  itself. 
Similarly,  the  lag  2,  3,  4,  5,  and  6  pairings  are  in  the  third, 
fourth,  fifth,  sixth,  and  seventh  diagonals.  (The  seventh 
"diagonal"  is  the  single  corner  element  since  the  only  lag  6 
pairing  is  the  first  and  seventh.) 

There  are,  however,  some  discrepencies  from  the  assumption 
that  the  correlation  coefficient  is  a  function  of  the  lag  only. 

The  QPF  set  and  the  WPF  set  without  the  outliers  display  some 
inconsistencies.  The  lag  1  and  the  lag  3  diagonals  have  somewhat 
broader  ranges,  in  the  QPF  set,  than  would  be  expected  if  the 
correlation  were  the  same  down  these  diagonals.  Similarly,  the 
WPF  set  without  the  outliers  has  a  bigger  range  than  expected  in 
the  lag  3  diagonal.  Also,  this  set  shows  an  odd  pattern  of 
negative  correlations  in  the  first  two  rows  and  positive  correl¬ 
ations  elsewhere.  This  is  an  indication  that  something  different 
is  happening  to  fatigue  crack  growth  rates  in  the  first  two 
increments  studied  here. 
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More  specific  conclusions  about  serial  correlation  are 
difficult  because  of  the  small  sample  sizes  and  the  limited  number 
of  increments  involved.  However,  from  the  broad  range  of  correla¬ 
tions  seen  in  Tables  3  and  4  it  is  clear  that  some  sort  of  serial 
correlation  exists.  Before  fatigue  crack  lifetimes  can  be 

adequately  modeled  from  FCGR ,  this  issue  must  be  resolved. 

4 

Virkler  et  al.  attempted  to  predict  fatigue  crack  lifetimes 
assuming  independence  of  FCGR  in  adjacent  increments.  Their 
technique  involved  fitting  a  different,  independent  lognormal 
distribution  to  the  FCGR's  in  each  increment.  They  used  these 
distributions  to  generate  a  set  of  fatigue  crack  growth  rates, 
which  were  numerically  integrated  to  get  a  FCG  life.  Their 
simulations  predicted  the  mean  FCG  life  very  well  but  were  off 
considerably  in  the  standard  deviation  of  FCG  lives.  The  pre¬ 
dicted  standard  deviation  was  approximately  one  third  the  actual 
standard  deviation  for  FCG  lifetimes  to  31mm. 
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SECTION  5 

THE  RELATIONSHIP  OF  FCGR  TO  STRESS  INTENSITY  FACTOR 


A  Paris  power  law  formulation  was  used  to  evaluate  the 
effect  of  stress  intensity  factor  on  FCGR.  The  form  of  the 
power  law  used  is: 


da 

dF 


C(iW 


m 


(1) 


where  C  and  m  are  constants.  The  stress  intensity  factor  was 

estimated  using  the  radial-thru-the-thickness  crack  formulas 

7 

suggested  by  Grandt  : 

K  =  cr/ira  [ - 0,8733 -  +  0.6762]  (2) 

0.3245  +  (|) 


where  R  is  the  radius  of  the  hole,  and  a  is  the  crack  depth 
measured  from  the  edge  of  the  hole.  No  attempt  was  made  to 
account  for  the  fact  that  the  crack  geometry  was  semielliptical 
growing  away  from  the  center  of  the  hole  and  down  the  bore.  It 
is  felt  that  the  Equation  2  formula  would  differ  from  the  actual 
crack's  stress  intensity  factor  by  a  multiplicative  constant 
although  this  constant  might  be  a  function  of  crack  length  if  the 
crack  shape  changes  as  the  crack  grows.  The  value  of  the  constant 
is  expected  to  range  between  0.7  and  1.0.  The  stress  intensity 
factor  given  in  Equation  1  is  determined  by  multiplying  the 
maximum  stress  in  the  fighter  stress  history  by  the  stress  inten¬ 
sity  factor  coefficient  obtained  from  Equation  2,  i.e. 


K  ~  o 

max  max 


(f)  =  {34  ksi)  •  (~) 


C  3 ) 


Procedures  similar  to  the  above  have  been  sucessfully  utilized 

g_g 

previously  to  describe  variable  amplitude  crack  growth  behavior 

By  converting  the  power  law  into  a  linear  equation  through 
a  log  transformation,  the  validity  of  the  power  law  model  can  be 
tested  with  linear  regression  procedures.  Linear  least  squares 
procedures  were  used  to  determine  the  power  law  constants 
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presented  in  Table  5  for  both  QPP  and  WPF  data  sets.  A  comparison 
between  the  individual  FCGR  data  points  and  the  power  law  deter¬ 
mined  to  describe  them  is  presented  in  Figures  5  and  6  for  the 
QPF  and  WPF  data  sets,  respectively. 

The  parameters  in  Table  5  indicate  that  the  power  law  is 
a  good  fit  to  the  data.  The  high  correlation  coefficients  (r) 
of  0.92  for  the  QPF  set  and  0.96  for  the  WPF  set  without  out¬ 
liers  indicate  a  strong  relationship  between  FCGR  and  stress 
intensity  factor. 

The  F  test  is  used  to  test  the  fit  of  the  model  and  is 
based  on  the  ratio  of  two  independent  estimates  of  the  variability 
of  the  FCGR  data.  One  of  these  estimates  measures  the  variability 
of  FCGR  after  the  model  has  been  fit  and  the  other  estimate  assumes 
no  model.  If  the  F  ratio  is  near  1,  almost  no  FCGR  variability 
is  explained  by  the  model.  If  the  F  value  is  larger  than  1,  the 
model  has  explained  some  of  the  variability  of  FCGR.  The  F  ratio 
has  an  F  distribution  whenever  the  model  does  not  explain  any  of 
the  FCGR  variability.  The  probability  of  getting  F  values  as 
large  as  those  in  Table  5  (P  value)  would  be  less  than  one 
hundreth  of  one  percent  if  the  model  did  not  explain  some  of  the 
FCGR  variability.  Therefore,  the  power  law  does  explain  some  of 
the  FCGR  variability,  or  in  other  words  is  a  good  fit. 

In  the  previous  section,  it  was  noted  that  the  QPF  specimen 
results  showed  a  faster  FCGR  than  the  WPF  specimens.  The 

i 

differences  in  FCGR  could  be  due  to  differences  in  crack  shape 

that  would  lead  to  different  stress  intensity  factors.  Table  6 

contains  the  thru-the-thickness  stress  intensity  factors 

“6  —6 

associated  with  crack  growth  rates  of  1.6  x  10  ,  4  x  10  ,  and 

-5  , 

10  inches  per  flight  hour  for  the  QPF  data.  The  stress  inten¬ 
sity  factors  for  the  corresponding  set  of  FCGR  values  were 
obtained  for  the  WPF  data  set  (without  outliers)  and  these  are 
also  presented  in  Tables  6  along  with  the  percent  difference  in 
stress  intensity  factor  values.  The  average  percent  difference 
for  the  WPF  set  is  about  8  percent  with  a  range  from  6.4  to  9.5 
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Table  5.  Results  of  Regression  Analysis 
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Figure  5.  A  Comparison  of  the  Power  Law  Fit  and  the  FCGR  Data  for  the  QPF  Set. 
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Set  Without  Outliers. 


percent.  The  difference  in  stress  intensity  factors  is  not 
surprising  since  there  could  be  as  much  as  a  thirty  percent 
difference  in  stress  intensity  factor  due  to  extremes  in  crack 
shape.  What  is  surprising,  however,  is  the  observation  that  the 
differences  became  greater  as  the  crack  grew  larger.  As  cracks 
get  larger,  they  grow  away  from  the  various  types  of  initiation 
scratch  patterns  and  should  grow  with  a  shape  that  is  independent 
of  their  starting  condition.  Thus,  what  we  really  expected  to 
observe  was  the  opposite  of  the  trend  shown  in  Table  6 . 
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Table  6.  Relative  Difference  in  K  as  a  Function 

c  T1  r*  r*  n  IH3.X 
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SECTION  6 
CONCLUSIONS 


The  fatigue  crack  growth  (FCG )  behavior  of  the  two  small 
crack  data  sets  studied  herein  is  similar  to  the  behavior  of 
large  cracks.  The  fatigue  life  distributions  for  these  small 
crack  data  sets  exhibit  a  lognormal  distribution.  The  coeffi¬ 
cients  of  variation  of  fatigue  crack  growth  rate  (FCGR)  for  these 

small  crack  data  sets  are  within  the  bounds  seen  in  other  studies 
dealing  with  longer  cracks.  The  FCGR  behavior  was  shown  to  be 
related  to  the  stress  intensity  factor  using  a  power  law.  All 
of  the  above  defined  traits  have  been  observed  in  large  cracks. 

The  two  sets  of  data  studied  were  slightly  different.  The 
fatigue  cracks  in  the  drilled  and  reamed  holes  set  grew  at  a 
slightly  faster  rate  than  the  fatigue  cracks  in  the  Winslow 
drilled  hole  set.  This  could  have  resulted  from  differences 
in  crack  shape.  Using  the  power  law  in  reverse  showed  that  the 
percentage  difference  in  stress  intensity  factors  that  would  give 
rise  to  the  same  FCGR  in  both  data  sets  ranged  between  6  and  10 
percent.  Differences  in  crack  shape  could  give  rise  to  this 
much  error  in  the  calculations  of  stress  intensity  factors. 

Before  an  adequate  model  of  FCG  lives  based  on  FCGR  can  be 
developed,  the  serial  dependencies  of  FCGR  must  be  accounted  for. 
Serial  correlation  was  observed  in  both  sets  of  data  analyzed  here ; 
but,  there  were  some  unusual  features  exhibited  for  cracks  below 
0.0125  inch.  Due  to  the  limited  number  of  data  points  specific 
conclusions  about  the  nature  of  this  correlation  cannot  be  drawn. 
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